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Abstract 

The Seebeck and Soret coefficients of ionically stabilized suspension of maghemite nanopar¬ 
ticles in dimethyl sulfoxide are experimentally studied as a function of nanoparticle volume 
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fraction. In the presence of a temperature gradient, the charged colloidal nanoparticles experi¬ 
ence both thermal drift due to their interactions with the solvent molecules and electric forces 
proportional to the internal thermoelectric field. The resulting thermodiffusion of nanoparti¬ 
cles is observed through Forced Rayleigh scattering, while the thermoelectric field is accessed 
through voltage measurements in a thermocell. Both techniques provide independent estimates 
of nanoparticle’s entropy of transfer as high as 75 meV.K^^ Such a property may be used to 
improve the thermoelectric coefficients in liquid thermocells. 


1. INTRODUCTION 

The thermoelectric effect (Seebeck effect) is known to influence the thermodiffusion behavior of 
charged colloidal suspensions with contributions coming from both electrolytes and charged col¬ 
loidal particles themselvesUnder a thermal gradient VT, the thermal drift of ionic species, 
i, induces concentration gradients Vnijni = —UiVT (Soret effect) and an internal electric field 
E — (Seebeck effect) Both Soret (a,) and Seebeck {S^g) coefficients depend on the East¬ 

man entropies of transfer 5, ,®® which characterize the interaction of species i with the solvent.!^ 
The absolute value of Si generally increases with the ion size. Large Soret effects have indeed been 
reported experimentally in various colloidal suspensions such as silica particles, DNA molecules, 
polystyrene spheres and magnetic nanoparticles (ferrofluids)l^^^ reflecting the large Eastman en¬ 
tropy of transfer associated with their equally large physical size (in the nm-/im range). The 
corresponding Seebeck coefficient, on the other hand, has rarely been studied.® 

Here, we investigate one such charged colloidal suspension, namely, ionically stabilized fer- 
rofluids in dimethyl sulfoxide (DMSO). Ferrofluids were chosen not only because of the high Soret 
coefficients®® of nanoparticles (NPs) but also because of their magnetic nature which may of¬ 
fer an additional control parameter (magnetic field®) on the thermodiffusive and thermoelectric 
behaviour. We determine independently the Eastman entropy of transfer of NPs S through (i) the 
Soret effect using forced Rayleigh scattering measurements and (ii) the Seebeck effect using a 
thermocell (see experimental section below). The values of S deduced from the two experiments 
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agree quantitatively and are almost three orders of magnitude higher than that of typieal ions in 
electrolytes ® Furthermore, we show that the thermodiffusion of NPs has a sizable influence on the 
Seebeck coefficient, an effect that may be used in liquid thermoelectric applications.^ 

2. EXPERIMENTAL SECTION 

2.1. Samples 

We have used ferrofluids based on well-known maghemite y-FeiOs nanoparticles (average diam¬ 
eter d = 6.7 nm and log-normal size distribution 0.38, determined from room-temperature mag¬ 
netization measurements) dispersed in DMSO. Nanoparticles are chemically synthesized first in 
water (see^for methods). This dispersion gives a ferrofluid of positively charged NPs with nitrate 
counterions, which are then replaced by perchlorate ones following the method described in.^i^ At 
the end of the process, DMSO is added instead of water to obtain electrostatically stabilized disper¬ 
sion with positively charged NP surface. The concentration of free perchloric acid in the solution 
was kept constant at i^l2.5 mM. This value was determined from the conductivity measurement 
in the supernatant of DMSO ferrofluid obtained after ultracentrifugation (60000 rpm, lh30) which 
separates the NPs from the solvent. 

2.2. Thermodiffusion measurements 

The Forced Rayleigh scattering technique used to extract the Soret and the NPs diffusion coeffi¬ 
cients is well described in Ref.^^The heating light (Hg arc lamp, lOOHz modulation) creates the 
optical image of a grid in the sample. Owing to the optical absorption by the NPs, a temperature 
grating is induced in the sample. Then a NP concentration grating settles due to the Soret effect 
in a few seconds at the spatial scale of ~ 50 /im. Both gratings are detected by the diffraction of 
a weakly absorbing test laser beam. As the gratings of the temperature and the NPs concentration 
evolve on timescales differing by orders of magnitude, this technique enables the use of a “two- 
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timescale” model.^The Soret coefficient is deduced from the temporal modulation of scattered 
intensity at a constant spatial modulation of NPs concentration. The diffusion coefficient is deter¬ 
mined by the relaxation time measurement of the concentration grating after the heating source is 
switched off. 

2.3. Thermoelectric measurements 

The Seebeck coefficient measurements were performed in a homemade thermocell consisting of a 
vertical, cylindrical Teflon cell (14 mm high and 6 mm diameter) with two ends sealed by sapphire 
windows, similar to the setup described in RefsThe ferrocene/ferrocenium (Fc/FJ) redox 
couple (2/4 mM respectively) was added to the sample in order to permit the exchange of electrons 
between the electrodes and the ferrofluid. These chemicals were purchased from SigmaAldrich; 
ferrocene (F^, 98%) and ferrocenium-tetrafluoroborate (FcBF 4 , technical grade) and used as re¬ 
ceived. The sample preparation is performed in a glovebox under a nitrogen atmosphere. It should 
be noted that the co-existence of the redox couple and nanoparticles did not change the redox po¬ 
tential (cf. Supporting Information I) or cause aggregation of NPs. The experiments are carried 
out between 30 °C and 50 °C (mean temperature) with the temperature difference between the 
two electrodes ATg/g^ ~ 4.3 °C (10 °C difference between the cell extremities). The open circuit 
voltage is: AV — —AATgiect, where A shall be referred to as the “thermoelectric coefficient”, to be 
distinguished from the Seebeck coefficient S^e- The thermoelectric voltage was monitored over 2 
hours between each temperature change. The diffusion coefficient D (see inset in Figure [^, gives 
the nanoparticle diffusion time T = /^/ {AkD) of the order of two days for our thermocell. Thus 
the observed electromotive force and thermoelectric coefficient A correspond to those of the initial 
stated 
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3. RESULTS AND DISCUSSION 

3.1. Theoretical considerations 

To analyse our experimental results, here we eonsider a eolloidal solution eontaining a eoneentra- 
tion n of eharged particles with a structural diameter d and a structural charge Z^tr (here, particles 
are positively charged). They are stabilized in a monovalent electrolyte solution A+B^. Some 
anions are condensed within the first solvation layers of the particles, partly canceling Zstr thus 
leading to an effective charge Z « Z^^^.I^The remaining anions- whose Coulomb binding en¬ 
ergy is smaller than ^gT-are free. The concentrations of free anions and cations are ri- and 
respectively. The electroneutrality writes: 


Zn -(-n+ — = 0. 


( 1 ) 


Under a temperature gradient, the particle current 7/ corresponding to the charged species i is: 




Ji 


-Di 


Vtli -h Hi 


Si 

ksT 


^T-m 



( 2 ) 


where n, is the particle density. The first term corresponds to Pick’s diffusion with coefficient D,, 
the second term, proportional to the “Eastman entropy of transfer” represents the thermal 
drift, and the last term is the electric drift in the presence of a local field E. The dimensionless 
number = kgTjif /eDi is proportional to the ratio of the electrophoretic mobility nf to the 
diffusion coefficient D, .I^For small point-like ions the Einstein relation is valid and is simply the 
ionic charge number Zi- The particle currents corresponding to small non-interacting monovalent 
ions are: 


7± = -T>± 


Vn±-h ^ Vr =F 

KbT KbT 


(3) 


For colloidal particles ^ is of the same order of magnitude as-but not equal to-the effective charge 
Z^ At large volume fractions, the interaction between NPs needs to be considered. This can be 
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described in terms of isothermal osmotic compressibility^ 0 = Vnpn where Vnp = TtcP{6 is 
the nanoparticle volume {cf. Supporting information II). The (j) dependence of the parameters S 
and ^ in eq ?? appears as: 


S = Soxi(l)) and ^=^ox{(j)). 


(4) 


Choosing a hard sphere model with Carnahan-Starling equation of state,l^;t^(^) becomes: 




1 + 40eyy + ^^eff “ ^^eff + ^^eff 


(5) 


where ^eff = ^{dns/d)^ represents an effective volume fraction corresponding to hard-sphere 
diameter dns = d-\- 2 Xd, where Xd is the Debye length. 

When the stationary state ('^0 is reached, each of the three currents expressed in eqs ?? and ?? 
vanishes. Combining these equations with the electroneutrality condition, eq ??, we obtain:!^ 




1 
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ZnS + 

+ ^Zn 


VT 


yfVT. 


( 6 ) 


Substituting this expression in the nanoparticle current (eq ??), we obtain: 


Vn 

n 


— {S-^eyf)VT^-aVT, 
ksT 


(7) 


where 

a = {S-^ey^')/kBT ( 8 ) 

is the Soret coefficient. For uncharged particles only the first term is present. In a series of re¬ 
cent papers, Wiirger and coauthors*have emphasized the importance of the second term 
in charged colloidal suspensions. Recent experiments by Eslahian et al^ on the salinity (elec¬ 
trolytes) effects on the thermodiffusion of polystyrene sulfonate beads appear to confirm these 
theoretical claims. 


6 









3.2. Thermodiffusion and thermoelectric analysis 

The Soret coefficient a, determined at different NP volume fractions 0 are reported in Figure 
with an inset showing linear variation of the diffusion coefficient: D{^) Ri 7.03 x 10 - 12(1 + 
66.02(j)) m^.s-i. The variation in a was analyzed using eqs ?? and ??, with the approximation 
(to be justified later): S » i-e.: 

l + i^i(l>,ff) + l)Zn/2 



Here, we define n = njn^, where (H+ ions) is kept constant. The nanoparticle structural 
diameter is d = 6.7 nm, the Debye length is A/j = 2.1 nm in a solution of 12.5 mM HCIO 4 in 
DMSO, with dielectic constant £ = 48, at room temperature. ^0 ~ 25 was estimated from the 
measurement of the electrophoretic mobility of a NP suspension at ^ = 0.05% using the laser 
Doppler velocimetry technique (NanoZS Malvern GB). The remaining unknown parameters Z 
and Sq are determined through the fit (solid line in Figure of the experimental data by eq ??. 
We obtain: Z Ri 30 and Sq ~ 68 meV.K-i, or equivalently 5o/i^o ~ 2.7 meV.K-i, one order of 
magnitude higher than typical values corresponding to electrolytes,^^ which justifies our previous 
approximation: S » ^qS+^^qS-. 



Nanoparticle volume fraction {% vol.) 


Figure 1: Soret coefficient a, and diffusion coefficient D (inset) as a function of NP volume frac¬ 
tion. Note that the measurements at higher concentrations (~ 4%) were performed on another set 
of DMSO based ferrofluids in similar ionic conditions. 
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In thermoelectric measurements, a temperature gradient Vr is established across a previously 
isothermal, homogenous electrolyte. While the bulk distribution of different species is still uni¬ 
form, an internal electric field, settles immediately within the fluid, resulting from charge ac¬ 
cumulations at the cell boundaries.This is due to the response of ions to thermal forces /,■ = SiVT. 
In open circuit conditions the total electric current: Jeiec — ZeJ -\-eJ+ — eJ is zero. Substituting 
eq ?? and eq ?? for the particle currents and taking into account the initial condition Vn, =0 we 

obtain :S1 




( 10 ) 


where ti = OijoT is the Hittorf transport number of ionic species i, i.e the relative contributions of 
its conductivity a, to the total conductivity: Oj = E/O'i- 

In a thermocell, where a reversible redox reaction occurs at the electrodes, the difference 
of electrochemical potential between the hot and the cold electrodes at initial state is: A/1 = 
—= eA'"''^Ar, with® (see Supporting Information III) for more detail): 






imt 




( 11 ) 


The first term A^^ = sp^+ —spc represents the redox reaction entropy at electrodes (i.e. the differ- 



Figure 2: A"“^ measured as a function of NP volume fraction (j) at cell median temperature of 30 
(circles), 40 (squares) and 50 °C (triangles). The solid lines are fits to the Eq. (??). See text for 
more explanation. 
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ence of the partial molar entropies of Fj and Fc) whieh remains eonstant through out the measure¬ 
ments. The seeond term arises from the initial eleetrie field as deseribed in eq ??. 

The measured as a funetion of nanopartiele volume fraetion (^), at different temperatures, 
are shown in Figure The measurements are repeated at least 5 times at eaeh eoneentration and 
temperature. The data dispersion is less than 4%. The A'"" values were found negative, as it ean 
be expeeted from the negative redox reaetion entropy of the Fc/F^ eouple,!^ and varied between 
-1.1 and -1.6 mV.K^^ The absolute value of A'""^ deereases with inereasing (j) as well as with the 
mean eell temperature. 

At a fixed mean eell temperature and eonstant HCIO 4 and F^/FJ eoneentrations, the variations 
in A'"'' stems from the term tS{^)/{^{^)e) of nanopartieles. The nanopartieles’ eontribution to 
the eleetrieal eonduetivity is: o = ^{<j))Ze^nD{<j))/ksT is less than a few pereent of the total 
eonduetivity Cr at 0 < 1%. The (j) dependenee of Ot was thus negleeted. Combined together, 
eq ?? ean be rewritten as: 


= ^"..>(0) + (12) 

ksTYnp Ot 

At room temperature, T>(0) is taken from the linear relation observed through Foreed Rayleigh 
seattering (Figure inset). Or = 35.5 mS.m^^ was measured at room temperature. Sinee the 
dominant temperature dependenee of both D and Oj arises from a same quantity i.e.: the inverse 
frietion eoeffieient l/r}{T) of the solvent, we shall take, as a first approximation, D/oj indepen¬ 
dant of T. The (j) dependenee of S(^) is obtained from the hard sphere model (eq ??). With Z = 30 
and d = 6.7 nm, the experimental results (Figure]^ are fitted to eq ?? to deduee Sq. The results are 
eompared in Figure to the Sq value determined from the foreed Rayleigh seattering experiments 
at 23"C. In the explored temperature range Sq determined from the Soret and the Seebeek eoeffi- 
eients measurements are ~ 75 meV.K^^, i.e. three orders of magnitude higher than the Eastman 
entropy of transfer of small usual eleetrolyte ions (e.g., 0.12 meV.K^^ for sodium ions in watei®). 
This observation supports the idea that the Eastman entropy of transfer play a major role in both 
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Figure 3: Eastman transport entropy per NP in infinite dilution limit Sq as a function of temperature, 
estimated by Seebeck (circle, square and triangle) and Soret (diamond) coefficients based models. 
The large error bars are mainly due to the uncertainty in the average NP size (see text). 

thermoelectric and thermodiffusive phenomena in charged colloidal solutions. 

4. CONCLUSION 

In this work, we have measured the thermoelectric coefficient of an ionically stabilized ferrofluid 
as a function of nanoparticle volume fraction and compared the results to the corresponding Soret 
effect measurements. As expected, both coefficients depend on the concentration of charged 
nanoparticles and the values of Eastman entropy of transfer, 5o, determined from both experiments 
are found to be in fair quantitative agreement. Our results lend strong support to the existing theo¬ 
retical models describing charged colloidal solutions’ thermoelectric and thermodiffusive proper¬ 
ties that both depend on Sq. Following the same rationale, one can postulate that the sign and the 
magnitude of Seebeck and Soret coefficients must depend on several experimental parameters: e.g. 
the relative importance of the Eastman entropy of transfer betwen nanoparticles and surrounding 
ions, the concentration of all charged species in the solution, and the surface charge of colloidal 
particles. These extensive parameters can be tuned experimentally to control the thermoelectric co¬ 
efficient A of charged colloidal suspensions, offering a new perspective in future liquid thermocell 
research. 
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